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Abstract. Wc present 4 fim ISAAC imaging and spectroscopy of the extremely metal poor dwarf galaxy SBS 0335- 
052, aimed at a better understanding of the dust in this low-metallicity galaxy. The 4 ^rn emission turns out to be 
very compact, confined to the brightest pair of Super Star Clusters (SSCs). The Ks-L' color is extremely red, and 
the L' emission is consistent with the extrapolation of the ISO mid- infrared spectral energy distribution (SED). 
From hydrogen recombination lines and a fit to the near-/mid- infrared SED, we confirm a visual extinction of ^ 
15 mag. Our data suggest that the sites of the optical and infrared emission are distinct: the optical spectral lines 
come from an almost dustless region with a high star formation rate and a few thousand OB stars. This region 
lies along the line-of-sight to a very dusty central star cluster in which there are more than three times as many 
massive stars, completely hidden in the optical. Prom the extinction, we derive an upper limit for the dust mass 
of 10^ Mq which could be produced by recent supernovae. 

Key words. Galaxies: dwarf; Galaxies: ISM; Galaxies: starburst; Galaxies: star clusters; Galaxies: individual: 
SBS0335-052 



1. Introduction 

With an abundance of .^0/41, SBS 0335-052 is the 
lowest-metallicity galaxy in the Second Byurakan Survey 
(Markarian et al. 1983), and the second lowest known af- 
ter I Zw 18. Because of its low optical luminosity {Mb = 
— 16.7), its small size (3-4 kpc diameter), and its strong 
narrow Hll-region-like emission lines, it is classified as a 
Blue Compact Dwarf (BCD) galaxy (Thuan et al. 1997). 
SBS 0335-052 hosts an exceptionally powerful episode of 
star formation that, according to Thuan ct al. (1997), oc- 
curs mainly in six Super-Star Clusters (SSCs) not older 
than 25 Myr. ISO observations of SBS 0335-052 (Thuan 
et al., 1999 - hereafter TSM) have revealed copious mid- 
infrared emission (the spectrum peaks at 14 /xm) and a 
spectral energy distribution (SED) that is well fit from 7 
to 17 /im by a heavily absorbed modified blackbody. From 
their fit TSM deduce an optical extinction in the range 19- 
21 mag and speculate that most of the star formation in 
SBS 0335-052 may be optically obscured. 

To better probe the extinction in SBS 0335-052, Vanzi 
et al. (2000 - hereafter VHTI) obtained a high spatial res- 
olution image in the Ks band and a near-infrared (NIR) 
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spectrum. Comparing these observations with HST im- 
ages and optical ground-based spectroscopy, they found 
no clear evidence for high extinction in SBS 0335-052. 
However, NIR colors suggest that a fraction of the if- band 
flux is produced by dust. Indeed that fraction is consistent 
with the mid-infrared SED, and a fit of the combined from 
2 to 15 /zm spectrum gave a visual extinction Ay ~ 12 
magnitudes (VHTI). 

To determine the nature of the dust, unexpected in 
such a low-metallicity object, we have targeted SBS 0335- 
052 at 4 /zm to obtain, for the first time, spatially resolved 
imaging at a thermal wavelength. At the same time, the 
low extinction at 4 /xm enables us to probe the ionized gas 
via the Bra recombination line. 

2. Observations 

We obtained a long-wavelength (2.7 - 4.2 /xm) low- 
resolution (R=360) spectrum of SBS 0335-052 with 
ISAAC at the ESO-ANTU (UTl). The observations were 
acquired in three different occasions on October 9, 2000 
(34 minutes integration) and on January 3 and 5, 2001 (1 
hour integration each). We used a 1" wide slit to match 
the already available NIR and optical spectra and a posi- 
tion angle PA=145°. The spectrum was acquired with an 
elementary integration time of 0.104 sec, while chopping 
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Fig. 1. ISAAC spectrum of SBS 0335-052. The horizon- 
tal dashed line shows the positions of zero flux, mak- 
ing evident the rise of the continuum toward long wave- 
lengths. The region of poor atmospheric transmission 
around 3.3 /^m has been removed. 



with the secondary and switching (North-South) between 
two beams with the telescope. The telescope was dithered 
along the slit after each beam-switching cycle. HR1891, 
a B2.5V star, was observed in the same way, in order to 
eliminate the telluric absorption of the atmosphere. We 
used the Eclipse package developed by ESO to reduce the 
spectra, following the standard procedures. The 2h34m 
combined spectrum, extracted with a 1" aperture, is dis- 
played in Fig. H Since the dispersion is 14A/pix, but the 
slit is almost 7 pixels wide, the original spectrum has been 
smoothed to the effective resolution dictated by the slit. 
The spectrum was flux-calibrated with the photometry of 
our L-band image described below. 

We used the same instrument and telescope to also ac- 
quire an image in the band at 3.8 fim (AA = 0.58 fim). 
The pixel scale was 0.0709"/pixel and the total integration 
time 30 minutes, with an on-chip integration of 0.104 sec. 
Like the spectrum, the image was acquired while chopping 
with the secondary and nodding the telescope between the 
two beams. The telescope was also dithered randomly af- 
ter each beam-switching cycle. We used the IRAF package 
for the data reduction, and relied on the telescope offsets 
to align the dithered images. The seeing FWHM measured 
from a star in the negative image is 6.9 pixels, or 0.49". 
The L' image is shown in Fig. |[ 

The L'-band photometry was calibrated with the pho- 
tometric standard HD 22686, obtained in the same way as 
the target, and assuming an L' magnitude of 7.20 (Elias et 
al. 1982). The growth curve of the object in the L' image is 
shown in the upper panel of Fig. the total L' magnitude 
is 14.1 ±0.2, shown as a horizontal dotted line. 




^ Because of its central wavelength 7^ 3.5 fim, it is more ap- 
propriate to call this band L' , which we will do hereafter. 



Fig. 2. i'-band ISAAC/VLT image with superimposed 
contours of the Ks image in red (VHTI), and those of the 
HST V image in blue (Thuan et al. 1997) (N up, E left). 
The SSCs denoted by Thuan et al. (1997) are also labeled. 



2.1. Photometry and colors 

SBS 0335-052 consists, as noted in the Introduction, of 
several SSCs superimposed on an extended blue under- 
lying gaseous envelope. The two brightest SSC groups, 
SSC 1-f 2 and SSC 4+5 (using the notation of Thuan et 
al. 1997), are separated by approximately 1.4". The regis- 
tration of the L' image is made difficult by the lack of field 
stars in the field-of-view of the direct image. Nevertheless, 
because of the N-S nodding direction, a relatively bright 
star ^ 26" E and 52" S appears in the image of the nega- 
tive beam. We have therefore used this star in the Digital 
Sky Survey (DSS) image to register the L' frame, since 
the star is too far away to appear in either the HST or the 
Ks image (VHTI). The astrometry relative to this star 
indicates that the L' source corresponds to the brightest 
SSCs, SSC H-2, although with ~ 0.5-0.7" uncertainty be- 
cause of the large pixels (1-7") of the DSS. Superimposed 
on our L' image in Fig. || are the contours of the HST 
image (blue) and the high-resolution Ks image (red). 

To verify the association of the L' source with 
SSC 1-1-2, we have plotted in Fig. ||the L' and Ks growth 
curves centered on SSC H-2 and SSC 4-1-5, the pairs of 
SSCs visible in the shorter-wavelength images. Because 
there was no L' emission at the SSC 4-1-5 peak, we de- 
rived the L' photometry by fixing the position to the Ks 
SSC 4-1-5 peak. It is clear from Fig. ^ that the radial trend 
of the L' emission is consistent with the Ks SSC 1-1-2 po- 
sition, since the L' emission from the SSC 4-|-5 is initially 
linear with radius, which is what is expected for no signal 
at the nominal center. (The SSC 4-1-5 curve then turns 
over at roughly the separation of the two SSC pairs.) 

Figure § (lower panel) shows the Ks-L' color obtained 
from combining our L' photometry with that from the Ks 
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Fig. 3. L'-band growth curve and Ks-L' cumulative pro- 
file of SBS 0335-052. The top panel shows the data points, 
with respective errors, and the horizontal dotted line the 
total magnitude of 14.11; a solid line illustrates the growth 
curve for a point source on the final image. Also shown in 
the top panel is the Ks growth curve derived from the high- 
resolution Ks image in VHTT The bottom panel shows 
the Ks-L' color. For all curves, the SSC l-)-2 component 
is shown as a solid line, and SSC 4-|-5 as a dashed one. 



image in VHTL The Ks-L' color is extremely red, ranging 
from 2.8 at the center to roughly 2 for the totalQ. After 
correcting the color for our redder L' filter (see Bessell 
& Brett 1988), we find K - L colors of 2.1 to about 1.5 
(these colors correspond to an L filter centered at 3.5 /zm). 
Such a red color is highly unusual in (non-Seyfert) extra- 
galactic objects, and is redder than all but one of the HII 
galaxies studied by Glass & Moorwood (1985). Indeed, 
the exception, NGC 5253, has K — L colors similar to 
SBS 0335-052, and has been called the "youngest starburst 
known" by Rieke et al. (1988). Because it is also a low- 
metallicity dwarf galaxy, we will use NGC 5253 below as a 
"benchmark" for further comparisons. These colors have 
not been corrected for the ionized gas emission (see Section 
H). 



3. The dust in SBS 0335-052 

The dust reponsible for the L' emission appears to be 
very compact, confined to a region at most 1.2" in diam- 
eter (see Fig. ^). At a distance of 54.3 Mpc (Thuan et 
al. 1997), this corresponds to approximately 300 pc. (This 
estimate is consistent with the physical extent of the ISO- 
emitting dust hypothesized by TSM.) The FWHM of the 
4 /im emission is roughly equivalent to the seeing FWHM, 
namely 0.5" or 130 pc. As seen in Fig. |[ the L' curve 
is very slightly more extended than a point source, but 
is much more compact than the Ks flux, which reaches 



its asymptotic value at a diameter of roughly 1 kpc (4"). 
(SBS 0335-052 is actually more extended than this at very 
low K surface brightness levels, see VHTI). 

3.1. Extinction 

Brackett a in emission is clearly detected in our spec- 
trum: we measure a flux of 9.0 ± 1 10~^^ erg/s/cm^ in an 
aperture of 1 " x 1.5 ". We can compare this value with the 
fluxes measured for Br7 by VHTI and H/3 by Izotov et 
al. (1997) in the same aperture, and use the ratio to de- 
rive the visual extinction (assuming a foreground screen) . 
Using the intrinsic line ratios given in Osterbrock (1989) 
appropriate for SBS 0335-052 (Case B, 20000 K, 400- 
500/cm'^: Izotov et al. 1997), together with the extinc- 
tion curve from Cardelli et al. (1989), we obtain from the 
H/3/Bra ratio Ay = 1.45 ± 0.11 mag. The Br7/Bra ra- 
tio gives Av — 12.1 ± 1.8. Izotov et al (1997) measured 
Av = 0.55 mag from Ha/H^, and Vanzi et al. (2000) 
Ay — 0.73 from HP/Bij. There is therefore a clear ten- 
dency of the extinction to increase with wavelength. 

With different assumptions, recombination line ratios 
can be used to infer either the interstellar extinction curve, 
or the geometry of the obscuring dust. In the latter case, 
longer wavelength line ratios have the virtue of probing 
deeper into embedded regions, an important advantage 
for dusty starbursts. Following Calzetti et al. (1996), we 
have calculated the color excess E{B — V) derived from the 
hydrogen recombination lines. As before, we adopted the 
Cardelli et al. (1989) expression for the interstellar extinc- 
tion curve, and the intrinsic line ratios from Osterbrock 
(1989). The color excesses are shown in Fig. |4| together 
with two dust models, a foreground dust screen and a ho- 
mogenously mixed slab of gas and dust. It is clear from 
the figure that a foreground screen model is highly inap- 
propriate both for SBS 0335-052 (x) and for NGC 5253 
(open circle, taken from Beck et al. 1996). It is also evi- 
dent that the model of homogeneous mixed slab predicts 
an increase of Ay with wavelength, qualitatively similar 
to our observations. Nevertheless, while NGC 5253 ap- 
pears to be well described by homogeneously mixed gas 
and dust, the low Ha/H/3 ratio in SBS 0335-052 together 
with the extremely high infrared line ratios are not well 
reproduced by a single value of r in the mixed-medium 
model.0 

Because the extinction derived from the optical line 
ratios is small {Ay = 0.55 mag), a more probable model 
consists of a nearly dustless region, responsible for virtu- 
ally all of the optical emission, lying in front of a highly 
obscured central knot. We exclude the possibility that 
the interstellar extinction curve could be grossly incorrect, 
since we have also used the Seaton (1979) and Landini et 
al. (1984) curves and obtained similar results. An extinc- 
tion law with a higher R = Ay / E{B — V) coupled with 
a mixed medium model could fit the data. Such ratios 



Because these curves are cumulative, such a trend implies a 
radial blueward gradient; the 2/im emission is more extended 
than that in L' . 



A clumped medium (Natta & Panagia 1984) is even worse 
at reproducing the observed line ratios. 
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Fig. 4. The color excess E{B~V) derived from the hydro- 
gen recombination line ratios; the data from SBS 0335-052 
are shown as x. Two simple models are also illustrated: 
the foreground screen model (labeled "Screen" ) and a ho- 
mogeneous mixture of gas and dust (labeled "Mixed"). 
The open circle shows data for NGC 5253 (Beck et al. 
1996). 



are typical of denser environments such as dark molecular 
clouds (e.g., Kim et al. 1994), but the apphcabihty of such 
a law to SBS 0335-052 is not clear. Therefore, in Section 
4 we will characterize the starburst in SBS 0335-052 using 
the simple two-component model outlined above, namely 
a low-Ay region in front of a region highly obscured. 

3.2. Emission 

We have added new 4 ^m and 2 /im points (VHTI) to 
the ISO spectrum of SBS 0335-052 (TSM), and fit the 
resulting SED. First, though, we had to estimate what 
fraction of the K and L' flux is due to dust. It turns 
out that the gas fractions determined in VHTI and be- 
low, together with reasonable estimates of stellar colors. 



constrain the relative contributions of ionized gas, dust, 
and stars from 1.6 to 4/im. The integrated Ks magnitude 
was used by VHTI to compare with the ISO observations 
obtained with a beam size larger than the galaxy; here, 
because of the compact 4 fim morphology, we adopt the 
VHTI Ks photometry for SSC 1+2 only. With the emis- 
sion coefficients given in Joy & Lester (1988), and our Bra 
observation, we find that 27% of the L' emission derives 
from the ionized gas continuum. Then assuming: i) no 
dust emission in H (VHTI); ii) stellar H - K ~ 0.0-0.2; 
in) stellar K — L' ^ 0.0 — 0.5, we derive a stellar fraction 
of ~ 37% at 2 /im, compared with 6% at 4 /in;j^. For sim- 
plicity, we have not considered any obscuration of the stel- 
lar component. With a 50% gas fraction at 2/im (VHTI), 
we therefore obtain a dust-emitting fraction of 13% at 
2 /im, and 67% at 4 /im, corresponding to 0.014 mJy at K 
and 0.37 mJy at L' . 

The combined SED from 2 to 17/im was fit with two 
modified blackbodies (MBBs); the cooler one is obscured 
by dust in a foreground screen with an extinction curve 
given by Lutz (1999). As shown by TSM, this curve pro- 
vides a much better fit than previous infrared extinc- 
tion curves, because of its significantly higher extinction 
in the 3-8 /im region. The emissivity was fixed to A^^'^ 
(TSM). Figure ^ shows the resulting fit; the warmer (un- 
obscured, shown by a dashed line) MBB has a tempera- 
ture of 459 ±20 K, and the cooler one (dot-dashed line) 
192 ±4 K, obscured by dust with Ay = 16.3 ± 0.5 mag. 
These results diff'er slightly from VHTI, but should be 
more reliable because of the additional constraint of the 
4 jim data point, and the use of a more refined extinction 
curve; also, in VHTI the emissivity was left as a free pa- 
rameter, and the 2 /im flux was higher because it included 
emission from the whole galaxy not just SSC 1-1-2 (see 
above). In any case, the flt is very similar to that found 
by TSM and confirms the high extinction found by them. 

Although our data suggest that the geometry is prob- 
ably more a homogeneously mixed medium than a fore- 
ground screen (see Fig. ^), the comparable values of the 
extinction {Ax/ Ay = 0.06 - 0.08 mag) in the MIR spec- 
tral region (Lutz 1999) make the shapes of the two fits 
very similar, but with Ay 4 to 5 times larger in the 
mixed medium model. We therefore adopted the fore- 
ground screen model since it gave fits of similar, if not 
better, quality than the mixed medium model. 

It was not possible to fit the SED with two MBBs 
absorbed ( "extincted" ) by the same amount, because the 
resulting fit falls off too steeply toward short wavelengths. 
Since the temperature and amplitude of the warm MBB 
are basically governed by only the 2 and 4 /im data points, 
an additional quantity such as extinction would not be 
constrained. However, the data do constrain the extinction 
for the warm MBB to be much smaller than that for the 
cool one (assuming the extinction curve is correct). Such 



The quoted values are the mean of the color range cited 
above. 
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Fig. 5. Near- and mid-infrared SED of SBS 0335-052 with 
residuals from the best fit. In the upper panel, the fit is 
shown by a solid line, while the dashed line shows the 
warm modified blackbody (449 K) which is not absorbed. 
The dot-dashed line shows the cool modified blackbody 
(193 K), and affected by a foreground screen of Ay = 16.3 
mag. The 2 /j,m point lying above the fit ( x ) illustrates 
the K value used in VHTI, roughly 10 times higher than 
that adopted here. The lower panel illustrates the frac- 
tional residuals which are largest around 12-14 /im, exactly 
where the extinction curve is most uncertain. 



a result is also consistent with the two-component model 
of the previous section. 

3.3. Mass 

The mid-infrared spectrum of SBS 0335-052 is unusual 
compared to other star-forming galaxies (e.g., Dale et al. 
2001), because of its positive 3-4 /im slope, its strong con- 
tinuum, and lack of spectral features (e.g.. Unidentified 
Infrared Bands: UIBs; Polycyclic Aromatic Hydrocarbons: 
PAHs) . Its dust properties are therefore not well described 
by standard dust models (Desert et al. 1990), especially 
regarding the very-small-grain component, and as a re- 
sult we cannot plausibly estimate the mass of the emit- 
ting dust. Nevertheless, the extinction is likely governed 
by the cooler dust, and we can estimate the mass of the 
obscuring dust from the Ay derived from our fit of the 
SED. 

With Av — 16.3 mag, and assuming N(III)/Av' ^ 
2 X IG^^cm^/mag and a gas-to-dust mass ratio of ~200, 
we derive a dust surface density of 1.2 A/0/pc^. This cal- 
culation assumes solar metallicity, but the metallicity cor- 
rection factors cancel out in the end. The derived dust 
surface density is similar to that of the model of Draine & 
Lee (1984) which, with a grain opacity of 3000 cm^/gm, 
is 1.6 Mq/pc^. If the extinction arises in a volume simi- 
lar in size to the emission region (diameter 1.2"), we ob- 
tain a dust mass of 10^ Mq. Alternatively, we can use 



the formalism of Spitzer (1978) to relate the amount of 
matter required to produce the observed extinction to the 
mean extinction curve, the density of solid material within 
the grain, and the grain dielectric constant (Aannestad 
& Purcell 1973). From the observed extinction, and as- 
suming the grain parameters of Draine & Lee (1984), 
we can therefore derive the mean volume density of the 
dust pd = 0.0065 Mq/pc"^, which is more than a factor 
of 20 higher than the Galactic value of 0.0003Mq/pc^. 
In a spherical region of diameter 1.2" as above, the dust 
mass becomes ~ 10^ Mq, consistent with the previous es- 
timate. These estimates are both upper limits since we 
have assumed a uniform distribution of the dust within 
the 1.2" region; the actual value may be less if the dust 
does not extend over the entire region. These upper limits 
are comparable with those obtained by TSM. 

4. Properties of the starburst 

Our measurement of the Bra line shows that a consid- 
erable fraction of the star formation in SBS 0335-052 is 
optically obscured, invisible even at 2 /im. Assuming the 
two-component geometry described in Section 3.1, we can 



adopt the (extinction corrected) II/3 flux measured by 
Izotov et al. (1997) to estimate the Br7 and Bra flux orig- 
inating in the same (low Ay) region. This exercise yields 
values that are respectively about 50 and 25% of what is 
observed. That is to say, 50% of the observed Br7 flux 
and 25% of the Bra comes from the (virtually) dustless 
foreground region which is the origin of the II/3 emission. 
The remaining 50% of the Br7 flux and 75% of Bra arise 
in the embedded dusty star cluster; in the optical we are 
only observing about 1/4 of the total ionized gas emission. 

The fraction of optically hidden emission can be used 
to estimate the extinction in the embedded cluster. We 
obtain Ay ^ 15, close to that inferred from the fit of 
the mid-infrared SED. The bright Brackett a line, cor- 
rected for Ay ^ 15, indicates that the total star forma- 
tion rate in the embedded cluster, instead of the optically- 
derived value of O.4M0/yr (Thuan et al. 1997), is more 
like ~1.7M0/yr. This last value is more than 40 times 
higher than that in the lowest-metallicity BCD I Zw 18 
(O.O4M0/yr). 

The number of massive stars in the embedded star 
cluster can be inferred from the the Bra luminosity. 
With 75% of the observed Bra flux coming from the 
obscured region, and correcting as above for Ay ~ 15, 
Lsra (embedded) = 4.2 x 10'^^ erg/s. Then, using the pre- 
scription of Guseva et al. (2000), we obtain approximately 
^ 14200 07 stars. There are more than three times 
as many massive stars in the embedded cluster as in the 
unobscured region ('^ 4040, with the observed Lnp lumi- 
nosity of 1.9 X 10"*° erg/s in the same aperture: VHTI). 
The total mass of the embedded cluster can be estimated 
by assuming an Initial Mass Function (IMF) ; adopting the 
IMF given by Scab (1998) gives a mass of 1.2x10^ M0 in 
the embedded cluster. Alternatively, with a Salpeter IMF, 
and lower and upper mass cutoffs of respectively 0.8 and 
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120 Mq (Schaerer & Vacca 1998), we obtain an embedded 
stellar mass of 6.6 x 10® Mq. The stellar density in the em- 
bedded cluster turns out to be ~ 85-150 Mq/pc^, which 
is relatively low compared to the SSCs in NGC 5253 with 
a stellar surface density of ^ 10^ Mq/pc^ (Calzetti et al. 
1997). 

The massive star cluster in the obscured central knot 
would be expected to host a significant number of Type II 
supernovae (SNe). Our Bra measurement, together with 
the lowest-metallicity models of Leitherer et al. (1999), 
imphes a SN rate of 0.004-0. GG6/yr, according to the stel- 
lar mass adopted (see above). We can therefore calculate 
the number of SNe expected to reside in the region of 
the SSCs, which turns out to be - 8000-12000, when in- 
tegrated over the lifetime of the burst after the onset of 
SNe. If each metal-poor SN produces, on average, 1 Mq 
of dust (Todini & Ferrara 2001), we would expect a dust 
mass on the order of 1O^M0, a factor of 10 lower than the 
upper limit derived from the extinction. This may suggest 
that the dust is not spread out uniformly over the 1.2" 
region. Also, the dust mass we infer from the extinction 
is highly uncertain because of the extremely low metallic- 
ity of SBS 0335-052. The physical conditions in the dense 
dusty medium of the embedded star cluster do not resem- 
ble those of the solar neighborhood, nor are they simi- 
lar to the cold dark clouds in the Galaxy. No dust mod- 
els exist for extremely low-metallicity environments; even 
the Magellanic Clouds modelled by Weingartner & Draine 
(2001) have metallicities about 4 times higher than that in 
SBS 0335-052. The grain properties such as the size distri- 
bution or the chemical composition may be radically dif- 
ferent in low-metallicity environments. It is also true that 
even in solar-metallicity contexts the canonical models are 
still a subject of debate. There is evidence, for example, 
that a more realistic size distribution may produce the 
same extinction with a lower dust mass (Kim et al. 1994). 
Also, some fraction of grains may be "fluffy" which would 
also increase the extinction per unit dust mass (Kriigel & 
Siebenmorgen 1994; Wolff et al. 1994; Mathis 1996; Snow 
& Witt 1996). 

SBS 0335-052, while unusual for such a metal-poor ob- 
ject, does have peers with similar properties, although at 
a higher metal abundance. NGC 5253 (~l/5 Zq) and 
Henize 2-10 {^1/10 Zq) are sub-solar metallicity low- 
luminosity dwarf galaxies (the former a dwarf elliptical, 
the latter a BCD) which host powerful central starbursts 
with SSCs (Calzetti et al. 1997; Conti & Vacca 1994), sim- 
ilar to SBS 0335-052. Both galaxies show high extinction 
from the Brackett recombination line ratios (Beck et al. 
1996, Kawara et al. 1989), but low extinction from the 
optical lines (see Fig. ^); they also show significant 10 /im 
absorption features (Kawara et al. 1989), and host signif- 
icant numbers of Wolf-Rayet stars in the central clusters 
(e.g., Guseva et al. 2000). Because of this and other ev- 
idence, the bursts are thought to be only a few Myr old 
(Beck et al. 1996; Calzetti et al. 1997; Beck et al. 1997; 
VHTI). In all three galaxies, the extinction is so high that 
optical measurements cannot be used to reliably study the 



embedded star clusters; even at 2 fim a significant fraction 
of the ionized gas emission is completely obscured. In the 
higher-metallicity objects (NGC 5253 and He 2-10), it ap- 
pears that the stellar clusters are born embedded in dust 
and then emerge after 2-3 Myr, becoming bright in the 
optical and the UV (Calzetti et al. 1997). Strictly speak- 
ing, such a scenario may not be appropriate for SBS 0335- 
052, since, if the present burst is the first one, dust could 
not have predated the current episode of star formation. 
The quantity of dust inferred from our measurements is 
roughly compatible with that from primordial SN mod- 
els, but the clusters themselves would have "polluted" the 
starburst. Either way, judging from SBS 0335-052 and its 
extremely low metal abundance, it is probable that much 
of primordial star formation could have occurred in deeply 
embedded dense star clusters. 

5. Conclusions 

1. We find direct evidence for a heavily absorbed central 
knot of star formation in SBS 0335-052. The extinction 
measured from the the observed Br7/Bra is Ay ~ 12 
mag, but line-ratio comparisons indicate that there is 
a virtually dust-free region along the line-of-sight to 
massive star clusters which are enshrouded in dust. 

2. The infrared SED has been fitted with two modified 
blackbodies, and the inferred extinction of 16.3 mag is 
roughly consistent with, but slightly larger than, that 
obtained from Br7/Bra. The total amount of obscur- 
ing dust is estimated to be around 10^ Mq. 

3. The star formation rate and number of massive stars in 
the embedded cluster in SBS 0335-052 are more than 
three times higher than inferred from optical measure- 
ments. Roughly 3/4 of the star formation in SBS 0335- 
052 occurs within a highly obscured embedded cluster. 

4. It is likely that the SNe produced by the current burst 
are responsible for the dust present in SBS 0335-052, 
although rough estimates of their dust production are 
lower than the upper limits for the dust mass inferred 
from extinction. 
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